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Abstract 

We present a nonperturbative dynamical study of e + e~ production in proton- 
nucleus and nucleus-nucleus collisions from AGS to SPS energies on the basis of 
the covariant transport approach HSD. For p + Be reactions the dilepton yield 
for invariant masses M < 1.4 GeV is found to be dominated by the decays of 
the rj, p, u) and $ mesons at all energies from 10 - 450 GeV. For nucleus-nucleus 
collisions, however, the dilepton yield shows an additional large contribution 
from 7r + 7r _ , K + K~ and np channels. Systematic studies are presented for the 
'free' meson mass scenario in comparison to a 'dropping' meson mass scenario 
at finite baryon density. We find that for 'dropping' meson masses the invariant 
dilepton mass range 0.35 GeV < M < 0.65 GeV is increased in comparison to 
the 'free' meson mass scenario and that the data of the CERES-collaboration for 
nucleus-nucleus collisions can be described much better within the 'dropping' 
mass scheme. We study in detail the contributions from the various dilepton 
channels as a function of the transverse momentum and rapidity of the lepton 
pair as well as a function of the charged particle multiplicity. Furthermore, 
various direct photon channels for S + Au at 200 GeV/u are computed and 
found to be well below the upper bounds measured by the WA80-collaboration. 



PACS: 25.75+r 14.60.-z 14.60.Cd 

Keywords: relativistic heavy-ion collisions, leptons 



*Work supported by BMBF and GSI Darmstadt. 



1 Introduction 



The question of chiral symmetry restoration at high baryon density is of fundamental 
interest since a couple of years [|I], [J, but a clear experimental evidence has not been 
achieved, so far. The enhancement of strangeness production as e.g. seen in the AGS 
data for the K + /tt + ratio [[| might be a signiture for such a transition ||, however, 
other hadronic scenarios can be cooked up to describe this phenomenon as well ||. 
Enhanced antikaon yields at SIS energies, as compared to transport studies using bare 
kaon masses || [7] , point in the same direction, but here the present knowledge on the 
elementary production cross sections close to threshold does not yet allow for a final 
identification. On the other hand, dileptons are particularly well suited for an inves- 
tigation of the violent phases of a high-energy heavy-ion collision because they can 
leave the reaction volume essentially undistorted by final-state interactions. Indeed, 
dileptons from heavy-ion collisions have been observed by the DLS collaboration at 
the BEVALAC [§, g, |T§ and by the CERES |TTJ, HELIOS 0, [lj, NA38 flTJ and 
NA50 collaborations [15 at SPS energies. 

Quite some years ago it has been found within microscopic transport studies at 
BEVALAC/SIS energies [[[6], [17) that above about 0.5 GeV of invariant mass (of the 
lepton pair) the dominant production channel is from 7r + 7r~ annihilation, such that 
the properties of the short lived p meson could be explored at high baryon density. 
The data available so far, however, did not allow for a closer distinction of the various 
models proposed. 

The recent data on e + e~ or p + p~ spectra at SPS energies, on the other hand, are 
more conclusive. The enhancement of the low mass dimuon yield in S + W compared 
to p + W collisions [|12| at 200 GeV/u has been first suggested by Koch et al. JTS[ 
to be due to 7r + 7r~ annihilation. Furthermore, Li et al. [HJ have proposed that the 
enhancement of the e + e~ yield in S + Au collision as observed by the CERES collabo- 



ration |Tl| should be due to an enhanced p-meson production (via it + it~ annihilation) 
and a dropping p-mass in the medium. In fact, their analysis - which was based on 
an expanding fireball scenario in chemical equilibrium - could be confirmed within 
the microscopic transport calculations in Ref. PD| . Meanwhile, various authors have 
substantiated the observation in Refs. [19, 20], that the spectral shape of the dilepton 
yield is incompatible with 'free' meson form factors ET|, [22], E3|. However, a more 
conventional approach including the change of the p-meson spectral function in the 
medium due to the coupling of the p, n, A and nucleon dynamics along the line of 



Refs. [^4[ [26j was found to be roughly compatible with the CERES data p(| ^7 
On the other hand, the dimuon data of the HELIOS-3 collaboration il2|l could 



too. 



only be described satisfactorily when including 'dropping' meson masses 0, |29] 



Though meanwhile there are a couple of hints on 'dropping' meson masses with 
baryon density and thus on a partial restoration of chiral symmetry for very hot and 
dense nuclear matter, a more systematic study of medium effects on the differential 
dilepton spectra - also with respect to the transverse momentum and rapidity y - 
appears necessary to optimize the experimental setups. 

Our paper is organized as follows: we start with a brief reminder of the covari- 
ant transport approach employed in our analysis in Sect. 2 and present the actual 
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expressions and form factors used for the evaluation of the differential dilepton mul- 
tiplicities. In Sect. 3 we investigate the question, at which bombarding energies one 
might find optimal conditions for high baryon densities at sufficiently large timescales 
in Pb + Pb collisions. Section 4 contains the comparison of our calculated dilepton 
spectra for p + Be and Pb + Au reactions with experimental data at SPS energies. 
Detailed predictions for dilepton production form AGS to SPS energies will be pre- 
sented in Sect. 5 for p + Be and Pb + Pb collisions within the 'free' and 'dropping' 
meson mass scenario. Section 6 concentrates on direct photon production in S + Au 
collisions at 200 GeV/u, while Section 7 concludes our study with a summary and 
discussion of open problems. 



2 The covariant transport approach 

In this paper we perform our analysis along the line of the HSDQ approach |3| which 
is based on a coupled set of covariant transport equations for the phase-space distri- 
butions fh(x,p) of hadron h [f|, |30|, i.e. 

{(n„ - n^ttf - Mtdpl) % + (n^x + M* h o;u s h ) d»} f h (x, P ) 

= Yl J d2d3d4...pG] 12 ^ 3i .J 4 (Tl + Tl 2 -Tl 3 -Tl 4 ...) 

hzhshi... 

x {fh 3 (x,p 3 )f h4 (x,p^f h (x,p)f h2 (x,p 2 ) 

- fh(x,p)f h2 (x,p 2 )fh ;i (x,P3)fh 4 ( x ,Pd} ■■■ ■ (!) 

In Eq. (H) Ufr(x,p) and U%(x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while 12^34...^ (II + n 2 — II3 — II4 . . .) is the 'transition 

rate' for the process 1 + 2— > 3+4+... which is taken to be on-shell in the semi classical 
limit adopted. The hadron quasi-particle properties in ([I]) are defined via the mass- 



shell constraint 30 



8{U»W-M?) , (2) 
with effective masses and momenta given by 

M* h (x,p) = M h + U*(x,p) 

n»(x,p) = p»-U%(x,p) , (3) 
while the phase-space factors 

f h (x,p) = l±f h (x,p) (4) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depend- 
ing on the type of hadron in the final/initial channel. The dots in Eq. ([!]) stand for 
further contributions to the collision term with more than two hadrons in the fi- 
nal/initial channels. The transport approach (HD is fully specified by U^(x,p) and 
U£(x,p) (/i = 0, 1,2,3), which determine the mean-field propagation of the hadrons, 
and by the transition rates G^G£ 4 (. . .) in the collision term, that describes the scat- 
tering and hadron production/absorption rates. 

1 Hadron String Dynamics 
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2.1 Hadron selfenergies 

The scalar and vector mean fields and U% for baryons are taken from Ref. |4[] and 
don't have to be specified here again, since variations in the baryon selfenergies within 
the constraints provided by experimental data were found to have no sizeable effect 
on the dilepton differential spectra. In the present approach we propagate explicitly 
- apart from the baryons (cf. 0) - pions, kaons, 77's, r]"s the 1~ vector mesons p, u, $ 
and K*'s as well as the axial vector meson a%. We assume that the pions as Goldstone 
bosons do not change their properties in the medium; we also discard selfenergies for 
the 7] and ?/-mesons in the present calculation, since we did not find any appreciable 
selfenergy effects in comparison to the experimental spectra available in the energy 
regime of interest here. 

The kaon selfenergies are described as in Ref. [Q following Kaplan and Nelson [ |3Tf , 
i.e. the kaon masses (K + , K°) are assumed not to be changed in the medium due to 
an approximate cancellation of attractive scalar and repulsive vector interactions, 
whereas the antikaons drop in mass according to 



m * R = m \ 1 - \ R — > (m u + m s ) w 0.16 GeV (5) 




with \x fa 0.16 as in Ref. [Q. In (|5|) m\ and p fa 0.16 fm~ 3 denote the vacuum 
antikaon mass and the saturation density of nuclear matter, respectively. Due to a 
lack of further knowledge about the K* vector mesons we assume the same scaling 
with baryon density as for the kaonsQ. 

The in-medium properties of the vector mesons p, uj and $ are modelled according 
to the QCD sum rule analysis by Hatsuda and Lee B3| as in our previous studies, i.e. 




m v = m° v 1 - Ay — > (m u + m d ) w 0.016 GeV (6) 



with Ay ~ 0.18 for p and to, while the $ meson is expected to scale with Ay ~ 0.025. 
In line with (|6]) the in-medium $ mass is limited to twice the strange quark mass, i.e. 
m$ > 2m s ^ 0.3 GeV. 

In addition to our previous works j|, |20], |28|, |33|] we now also include the axial 



vector meson a\, i.e. the chiral partner of the p meson. Since the scaling of the a± 
mass with baryon density is quite a matter of debate and we only know that m* ai = m* 
in the chiral limit, we assume the a\ to scale in mass in the same way as the p meson 
according to Eq. (|6|) in order not to introduce any new parameter. 

In the following, the 'free' meson mass scenario will employ the limit A^ = Ay = 
A ai = 0, whereas the 'dropping' mass scenario is described by Eqs. (|5]) and @, re- 
spectively. 



2.2 Dilepton channels 

In this analysis we calculate dilepton production by taking into account the contribu- 
tions from the Dalitz-decays rj — > / ye + e~ , u — > 7r°e + e~, 77' — > 7e + e~, a\ — > 7re + e~and 

2 This working hypothesis needs to be controlled in future. 
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the direct dilepton decays of the vector mesons p,u>, $ where the p, $ and a\ mesons 
may as well be produced in inr, KK and np collisions, respectively. 

In case of a perturbative treatment^ of the channel 7t + tt~ — > p° — > e + e~ the cross 
section is parametrized as in Refs. [|16|, |35|, |36| by 



7T~ r 7r — >e" r e 



47r a 2 / 4m?. . l0 



3 M 2 V M 2 
where the 'free' form factor of the pion is approximated by 

4 

' F ^ M ^' = (M 2 — 77i p 2 ) 2 + m 2 r 2 ' 
In Eq. (0) M is the dilepton invariant mass, a is the fine structure constant, and 

m p = 775 MeV, T p = 118 MeV. 



The cross section for K + K annihilation is parametrized as [07 



a K+K -^ e+e -(M) = -j (-) f " l^(M)| 2 , (9) 

where the 'free' form factor of the kaon is approximated by 

1 777 

with 

m$ = 1020 MeV, r$ = 4.43 MeV. 

The cross section for dilepton production in n + p~ — ► $ — ► e + e~, 7r~p + — ► $ — > 
e + e~ scattering can be represented as 

1 B 

&np^<l>~,e+e-(M) = - ^ G K + K -^^ e + e - (11) 



3 B K + K 



with = 0.13, B K+K - = 0.49. 
The 7] Dalitz decay is given by 



dT^e+e- _ 4a r^ 27 4m 2 y /2 ^ | ^ m\ N 



rfM 3tt M V M 2 ; V M 2 / 

/ M 2 \ 3 

x(^l-— J |F^ 7e+e -(M)| 2 , (12) 
where the form factor is parametrized in the pole approximation as 

^M 1 -^) ■ (l3) 

3 The perturbative treatment is used to test the dynamical scheme described below in case of the 
bare p-mass scenario, only. 
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with the cut-off parameter A„ ~ 0.72 GeV. 



Similarly, the u Dalitz-decay is |38 
dT 



2a T a; ^ 7r o 7 ( 



Am 



2\ 1/2 




' M 2 



(M)p 



(14) 



where the form factor squared is parametrized as 



.(M)p 



A 4 , 



(15) 



with 



A w = 0.65 GeV, = 75 MeV. 



For ?]' — > / ye + e~ we use a similar expression as Eq. ([12]) . However, in this case the 
pole approximation (|T3|) is no longer valid since the vector meson pole occurs in the 
physical region of the dilepton spectrum (M < m^). Instead we use a form factor of 
the form (|l^) with 



0.75 GeV, 



0.14 GeV, 



which reproduces the experimental data from |38|] reasonably well. 
The direct decays of the vector mesons u, $ to e + e~ are taken as 



da 
dM 2 



(M) 



m v Ti 



7r(M 2 



m v ) 2 + Ul^/Tlr 



(16) 



T 1 V 



with T^ e+e -/T^ = 2.5 x 10~ 4 and T^+^/T^ = 7.1 x 10~ 5 . 

The dilepton channels rj — > e + e~, — > 7r°e + e~, — > e + e _ , $ — > e + e~, 7rp — ► 
$ — > e + e~ and K + K~ —>■ $ — ► e + e~ are treated perturbatively and are computed at 
the end of the transport calculation due to the 'long' lifetime of the mesons 77, 77', a;, $. 
The decays of the short-lived p and a\ mesons, however, have to be treated explicitly 
since these mesons change their properties rapidly (in case of the dropping mass 
scenario) and decay during the expansion phase of the system. 

The mesons (p, a{) stemming from a string decay with invariant mass m v at baryon 
density ps according to Eq. are selected by Monte Carlo according to the Breit- 
Wigner distribution: 



f(M) = N v 



Mm v T* v 



vr (M 2 



rn 



*2\2 



+ m*JT 



(17) 



v v 



where Ny guarantees normalization to unity, i.e. / f(M)dM = 1. The width T V (M) 
is determined from 



r 'vM=Hw) t> (18) 

where T is the full width at the mean resonance energy, q and qy are the pion 
three-momenta in the restframe of the resonance with mass M and m v , respectively. 



6 



Apart from string decay the mesons (p, ai) are abundantly also created from mi or 
up collisions, respectively. For the a\ formation cross section in the reaction 7r + p~ — ► 
a 1 ,7r _ p + — > a\ or for the p° cross section the reaction 7r + 7r~ — > p c 
Wigner form |3{], 40]: 



we use the Breit- 



2J V + 1 4vr 



2S + 1 k 2 (s 



*2p*2 



m 



*2\2 
V I 



+ m* 2 Y* 2 ' 



(19) 



where k is the pion momentum in the center-of-mass of the produced meson V = 
(p, ai), s is the invariant energy squared while Jy stands for the spin of the produced 
meson and 5* for the spin of the collision partner of the pion, respectively. In the time 
reversed processes the vector mesons of actual mass M may decay in each timestep 
according to the probability 



P = exv(-r* v (M)/jAt) 



(20) 



where At is the actual timestep size and 7 the Lorentz factor of the resonance with 
respect to the calculational frame. 

The p° decay to e + e~ with invariant mass M is calculated by integrating the 
equation (using the mass bin AM) 



AMdt 



1 N p0 (M,t) 



E 



events i=l 



(M) 



1 



AM 



(21) 



in time with 

I>^ e+e -(M) = 8.8 x KT 6 M, (22) 

where iV p o (M, t) is the number of p° mesons of mass M at time t in the calculation. 
The factor 8.8 x 10~ 6 stems from the measured width of the p° to e + e~. It can be 
shown |J7| that the method described above leads to the same result as Eq. (|7]) if the 
p° meson does not change its properties in time. We note that by treating explicitely 
the p formation by tttt collisions the perturbative channel (|7|) has been switched off 
to avoid double counting. 

The ai decay to ire + e~ with invariant mass M is calculated in analogy by inte- 
grating the equation 



e^e 

dMdt 



E 

events i=l 



dT ai . 



dM 



(rn* ai ,M) 



(23) 



in time with 
dT 



a\— >7re~ l "e" 



dM 



2a T 



3tt M 



M 2 



Am 2 e 
2 



m 



*2 
"1 



mt 



1/2 



1 + 2 



M 2 



AmfM 2 

u, l 



m; 



*2 



m 



2\2 



3/2 



(24) 
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with r ai ^ 7r7 = 6.4 x 10~ 4 GeV, while N ai (t) is the number of a% mesons at time t in 
the calculation. 

We note that we discard baryon-baryon (BB), meson-baryon (mB) and meson- 
meson (mm) bremsstrahlung channels as well as the Dalitz decays of the baryon 
resonances since their contribution was found to be small in Refs. |f20| , f28f . 

Before going over to the actual calculations for e + e~ spectra, we start with some 
more general analysis of central heavy- ion collisions from 1 - 200 GeV/u. 



3 Optimizing for high baryon density 

In order to probe the restoration of chiral symmetry at high baryon density in nucleus- 
nucleus collisions, one has to perform experiments with heavy nuclei (e.g. Pb + Pb) 
and optimize the beam energy to achieve a large volume of high baryon density for 
a sufficiently long time. In this respect central collisions of Pb + Pb have been 
investigated within the transport approach specified above and the 'stopped' baryon 
density p s B (t) - including only baryons with rapidity \y\ < 0.7 in the cms - has been 
computed in a central cylinder of the volume V = 7rR 2 A z / / -f cm with A z = R = 4 fm, 
while 7 cm is the Lorentz factor in the nucleus-nucleus center-of-mass system. Since 
we are interested in high baryon densities above some value p m i n for long times, we 
consider the quantity 

dt (p%(t) - Pmin ) Q(p s B (t) - p min ), (25) 

which should serve as a useful guide in the optimization problem. The quantity F 
(Eq. (|25|) ) is displayed in Fig. [I] for central collisions of Pb + Pb from 1 - 200 GeV/u 
for different values of p m i n from 2p — 5po (po ~ 0.168 fm~ 3 ). Accordingly, optimal 
bombarding energies for baryon densities above 4p should be around 20 - 30 GeV/u 
in order to explore the properties of an intermediate phase, where the chiral symmetry 
might approximately be restored and the hadron masses (except for the Goldstone 
bosons) might be close to their current quark masses m q + m 5 . However, also lower 
bombarding energies (2 - 10 GeV/u) are seen to qualify for studies of a partial restora- 
tion of chiral symmetry since sizeable space-time volumes with baryon densities above 
2 - 3 po can be achieved. 

Whereas the HSD approach has been shown in Ref. Q to reasonably reproduce 
hadronic spectra and rapidity distributions for heavy-ion collisions at SIS and AGS 
energies, a definite proof for heavy systems at SPS energies is still lacking though 
quantitative predictions for the baryon rapidity distribution in Pb + Pb collisions at 
158 GeV/u have been given in Ref. [|[]. In order to demonstrate, that the analysis 
presented in Fig. |l] is meaningful also at higher bombarding energy, we compare in 
Fig. ^] the preliminary rapidity distribution of negative hadrons (essentially tt~ , K~ 
and p) from NA49 |[| for central Pb + Pb collisions at 160 GeV/u with the HSD 



results. Here the solid line corresponds to a calculation (at b = 2 fm) including 
the 'dropping' meson masses from Eqs. @, (H), whereas the dotted line results from 
a calculation with bare meson masses. The broadening of the rapidity distribution 
around midrapidity (y w 3) in the 'dropping' mass scenario is due to pions from p 
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and u decays, which are produced with a wider distribution in rapidity in this case. 
The quantitative agreement with the data, however, indicates that both scenarios - 
i.e. 'free' and 'dropping' meson masses as described above - are compatible with 
the present preliminary data. We note that the open squares in Fig. |2| have been 
obtained by reflecting the full squares (from NA49) at midrapidity. The fact, that the 
calculated spectrum is not fully symmetric with respect to midrapidity, is caused by 
the finite statistics of the computation, which was performed with 150 parallel runs 
in this case. 

Since the system Pb + Pb at 160 GeV/u is explored experimentally at the SPS 
in great detail, it is advantages to have a look at the space-time evolution of the 
baryon one-body density for central collisions. In this respect we show in Fig. [I] the 
space-time evolution of baryons for this system at b = fm: (l.h.s.) countor plot of 
the baryon density distribution in coordinate space pb(x,V = 0,z;t), (middle colum) 
contour plot of the baryon momentum distribution PB{Px,Py = 0,p z ;t), (r.h.s.) the 
phase-space distribution 

f{z,p z ;t) = (2tt)~ 2 / dr ±dp±_ f b {r±,z,p±,p z ;t), (26) 
b J 

where J2 denotes a sum over all baryon species. In the time evolution of the density 

b 

distribution p B (x,y = 0,z;t) we explicitly mention the short phase of high baryon 
density from about 5-8 fm/c as well as the sizeable fraction of 'spectators' from the 
nuclear corona. The time evolution in momentum space (middle colum) shows that 
the system reaches its final distribution within a few fm/c, however, is far from the 
kinetic equilibrium in the baryon degrees of freedom, which would be reflected by a 
spherical distribution here. It clearly indicates a dominant longitudinal expansion of 
the system, which is much more pronounced than at AGS energies (cf. Fig. 21 in H). 



4 Comparison with experimental dilepton data at 
SPS energies 

The covariant transport approach HSD has been applied already to the analysis of 



dilepton spectra for p + Be and Ca + Ca collisions at BEVAL AC/SIS energies [[12 
for p + Be, p + Au, and S + Au collisions at SPS energies |20| and to dimuon 
spectra for p + W and S + W at 200 GeV/u @ in Ref. gj|. Whereas the dilepton 



spectra for p + A reactions at high energy could be well described by the mesonic 
decays in line with the experimental analysis |L1[], the S + A data could only be 
satisfactorily reproduced within the 'dropping' mass scenario (cf. also [[B5], |2D[). Here, 
we include additional reaction channels (e.g. the rf, a\ degrees of freedom) and present 
a systematic study with respect to differential spectra in the transverse momentum of 
the lepton pair, the rapidity distribution as well as the dilepton yield as a function of 
the charged particle multiplicity at midrapidity. Before doing so, we present the actual 
status of our calculations in comparison with presently available data on dilepton 
production at SPS energies. 
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4.1 Differential dilepton spectra for p + Be and Pb + Au at 
SPS energies 

As an example for dilepton spectra at SPS energies Fig. |] shows the spectral decom- 
position as a function of the e + e~ invariant mass M for p + Be at 450 GeV/c in 



comparison to the data of the CERES collaboration |jlT| . The e + ,e~ acceptance cuts 
in pseudo-rapidity (2.1 < rj < 2.65), a cut of the transverse e + and e~ momenta for 
p T > 0.05 GeV/c as well as a cut on the opening angle of the e + e~ pair (0 > 35 mrad) 
are taken into account. Furthermore, the experimental mass resolution has been in- 
cluded in evaluating the final mass spectrum, which is normalized by the number 
of charged particles drich/drj in the pseudorapidity bin 2.1 < r] < 3.1. As can be 
extracted from Fig. ^, the spectrum for p + Be can be fully accounted for by the 
electromagnetic decays of the rj,r]',ai and vector mesons p°,u> and contributions 
from meson-meson channels (7r + 7r~, K + K~, np) are of minor importance here. We 
note that in our present analysis the a\ and rj' Dalitz decays have been taken into 
account in addition to Ref. p0| . The rf contribution is more sensitive to the electro- 
magnetic form factor then rj or uo Dalitz decays because the vector meson pole shows 
up in the region M < m v r. However, the contribution of the r( Dalitz decays is not 
essential compared to the other channels; the a\ contribution is practically negligible 
for p + Be. 

The situation changes quite dramatically when going over to nucleus-nucleus col- 
lisions. In Fig. |^ we compare the results of our calculation for the differential dilepton 
spectra for Pb + Au at 160 GeV/u at b = 5 fm with the preliminary experimental 



data |]43|j . Contrary to p + Be reactions, a cut of the transverse e , e~ momenta 
Pt > 0.175 GeV/c has been taken in line with the experimental acceptance cut. For 
Pb + Au at 160 GeV/u (and semicentral collisions) the dominant yield for invariant 
masses 0.3 GeV < M < 0.7 GeV stems from 7r + 7r~ annihilation (cf. Fig. |5|). Also 
in the $ mass regime about 1 GeV there is a large contribution from K + K~ and up 
annihilation to dileptons for both scenarios: with bare meson masses (upper part of 
Fig. [5p and with in- medium meson masses (lower part of Fig. [5]). Whereas most of 
the processes (Dalitz and direct decays) occur in the vacuum at zero baryon density, 
the 7T7T — > p° — > e + e~ and direct p° (from baryon-baryon and meson-baryon collisions) 
decay still occur at finite baryon density such that a dropping p mass also leads to 
a shift of the respective contribution to lower invariant masses M. In Fig. I both 



scenarios are compared to the preliminary data of the CERES collaboration |^3| ; due 
to the present statistics, however, there is no unique conclusion since the calculation 
with bare meson masses (upper part) also describes the data except for one point 
at 0.6 GeV (cf. also Ref. p3|). On the other hand, the present preliminary data 
match well with the calculation including the in-medium meson masses. We note 
that the comparison with the data in Fig. |5] has been performed for 6 = 5 fm, because 
for this impact parameter the charged particle multiplicity dn ch /drj ~ 260 as for the 
experimental normalization (see below). 
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4.2 Dilepton yield versus charged particle multiplicity 



The number of charged particles dn c h/dr] in the pseudorapidity bin 2.1 < rj < 3.1 
for Pb + Au at 160 GeV/u is shown in Fig. The open circles are the result 
of our computations with 'free' meson masses, while the solid circles correspond to 
calculations when including the in-medium modifications of the meson masses. In both 
cases the charged particle multiplicity decreases with impact parameter practically 
linearly. For peripheral collisions there is no essential difference between both schemes 
as expected. For central collisions, however, the charged particle multiplicity in the 
'dropping' meson mass scenario is slightly larger due to the reduction of the vector 
meson and antikaon production thresholds, which enhances the respective particle 
formation cross sections at high baryon density. Especially the subsequent decay of p 
and io mesons to pions leads to a larger number of pions in the final expansion phase. 
We note, that due to the conservation of energy and momentum in each production 
event the enhanced number of vector mesons and antikaons at finite baryon density 
goes along with a lower number of those mesons, that do not change their quasiparticle 
properties in the medium (ir,r),r]' etc.). 

Including the CERES acceptance cuts and mass resolution as described above, we 
show in Fig. |7] the dilepton yield integrated over the invariant mass range 0.3 < M < 
1.0 GeV, 

dN/dr] = 1,0 dM dn e+e ~/{dMdr]) 

1 ' 0.3 GeV 1 ' 

as a function of the charged particle multiplicity dn c h/dr] without (open circles) and 
with (solid circles) in-medium mass modification. At small charged particle multi- 
plicity, which corresponds to very peripheral collisions, the integrated dilepton yields 
coincide for both cases. With decreasing impact parameter the average baryon den- 
sity and especially the pion density increases; as a consequence the contribution from 
pion annihilation to p° and subsequent decay to dileptons becomes larger. Since we 
gate on dileptons above the 7r + 7r~ annihilation threshold, also the integrated dilepton 
spectra increase with dn c h/dr]. Using 'free' meson masses we reach some plateau for 
low impact parameter which implies that the irn annihilation contribution, divided by 
the pion density, becomes approximately constant. However, for the 'dropping' mass 
scenario the absolute dilepton yield above 0.3 GeV is smaller because the directly 
produced p-mesons (at high initial baryon density) also 'shine' in the invariant mass 
regime below 0.3 GeV (cf. lower part of Fig. |J). Furthermore, due to an initially higher 
vector meson density the initial pion density (due to energy conservation) is reduced 
as compared to the 'free' mass scenario and the corresponding pion annihilation con- 
tribution is also lowered to some extend. All effects together lead to an approximately 
linear increase of the integrated dilepton yield with the charged particle multiplicity 
in the 'dropping' mass picture. Experimental data with sufficient statistics should 
allow to disentangle the two schemes or disqualify the hadronic scenario as employed 
in the HSD transport approach. 
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5 Systematics of dilepton production from AGS to 
SPS energies 

In this section we present a systematic analysis of various dilepton observables - the 
differential dilepton spectra, rapidity and transverse momentum distributions, the 
average transverse momentum - for p + Be and Pb + Pb collisions from AGS to SPS 
energies for the 'free' and 'dropping' mass scenarios. 

5.1 Differential dilepton spectra 

We first examine the differential dilepton spectra - integrated over rapidity and trans- 
verse momentum - with respect to their 'cocktail' decomposition. The respective 
spectra for p + Be at 10, 50, 450 GeV/u and for Pb + Pb at 10, 50, 160 GeV/u for 
b = 2 fm are shown in Figs. |8] - [[(]. For all cases we include a mass resolution of 



AM = 10 MeV, which can be expected for future dilepton detector systems |[44|| . As 
seen from Figs. ^| - |H], there are no dramatic differences in the relative contribution 
of the various dilepton channels; the total yield (and especially the <3> decay) increase 
with bombarding energy quite smoothly. In case of central Pb + Pb collisions the 
dominance of the 7r + 7r~ annihilation component is most pronounced at 10 GeV/u in 
both scenarios and decreases with bombarding energy. 



Since experimentally only the total spectra can be observed, we show in Fig. |TT 
the sum of all contributions for the 'free' (dashed lines) and 'dropping' mass scenario 
(solid lines) for central collisions of Pb + Pb at different bombarding energies to 
demonstrate the influence of in-medium effects for the mesons. For a mass resolution 
AM = 10 MeV one can expect to observe not only the "usual" enhancement of the 
dilepton spectra by about a factor of 2 at invariant masses 0.3 < M < 0.6 GeV, but 
also a sharp drop of the spectrum above the u mass by a factor 4-6 due to the shift 
of the p contribution to lower invariant masses. Furthermore, the peak from the u 
meson becomes more pronounced since the 'background' from the p decay (either from 
ttB, BB or 7T7T collisions) is significantly reduced. In the $ mass region, furthermore, 
we find a small increase of the yield for the 'dropping' mass scheme at all energies 
from 10 - 160 GeV/u. Again, the relative modifications of the total spectrum are 
most pronounced at 10 GeV/u. 

5.2 Transverse momentum distributions 

In this subsection we explore if the transverse momentum distribution of the lepton 
pair might give some further criteria to distinguish between the different scenarios. 
In this respect the transverse momentum distributions - integrated over rapidity and 



the invariant mass 0.4 < M < 0.7 GeV - are shown in Figs. [L2| - [14| for p + Be at 10, 
50, 450 GeV/u and for Pb + Pb at 10, 50, 160 GeV/u for b = 2 fm using g ± -bins of 
50 MeV/c. 

The main contribution in the invariant mass region 0.4 < M < 0.7 GeV for p + Be 
comes from the rj, uj Dalitz decays according to Fig. ||. The same decomposition can 



be observed in Fig. ITI3, where we display only the dominant channels. Performing 
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an exponential fit to the sum of all contributions (thick solid lines in Fig. |12|) for 
M > 0.4 GeV the following slope parameters can be extracted from our calculations: 
125 MeV at 10 GeV/u, 160 MeV at 50 GeV/u, and 185 MeV at 450 GeV/u. 

Contrary to p + Be reactions, the dominant channels in the mass range 0.4 < 
M < 0.7 GeV for central Pb + Pb collisions are the pion annihilation, u Dalitz 
decay and direct decay of vector mesons. As shown before, including the in-medium 
'dropping' masses we obtain a shift of the p meson contribution (cf. Figs. |9|,PH) to 
smaller invariant masses and one might expect same enhancement of the transverse 
momentum distribution at lower q±_ in this case. This effect is demonstrated in Fig. [TB| 
where the transverse momentum spectra for the in-medium mass scheme (solid lines) 
are slightly softer than for the 'free' mass case (dashed lines). This tendency can 
also be seen from the slope parameters extracted for > 0.4 GeV/c: 145 MeV at 
10 GeV/u for the 'free' meson mass scheme - 130 MeV, when including the in-medium 
mass modification; at 50 GeV/u - 170 MeV and 155 MeV, respectively; at 160 GeV/u 
- 195 MeV and 180 MeV, respectively. 

In order to clarify the origin of these differences, we present the channel decompo- 
sition for central Pb + Pb collisions at 160 GeV/u in Fig. [14| for both scenarios. The 
pion annihilation and u Dalitz decay contribution, which are dominant in the invari- 
ant mass region considered, become slightly softer in the 'dropping' mass scenario. 
We note, however, that the transverse momentum spectra do not differ sizeably; ac- 
cording to the authors point of view they do not qualify very much for disentangling 
the different schemes. 



5.3 Average transverse momentum < qr > (M) 

The dilepton average transverse momentum as a function of the invariant mass, 
< qx > (M), was already studied experimentally for hadron-hadron collisions a cou- 
ple of years ago (cf. the reviews |4^, (461), where an increase of < qr > with M and 
bombarding energy was observed. 

In our analysis the average transverse momentum for the channel k is defined as 

Ew*(q Ti ,Mi) q Ti 
< q T > k (M) = ' — -, (28) 

k i 

where wf(qTi,M) is the probability of the individual dilepton event % (cf. Section 
2.2). The total < qx > is the sum of the individual channels: 

< q T > (M) = ]T < q T > k (M). (29) 



The results of our computations are shown in Fig. [15] for p + Be at 10, 50 and 450 GeV 
including the contributions from the individual channels. Due to mass bins of 50 MeV 
the shapes from the vector mesons u and $ are not very pronounced. However, one 
can see the main tendency: the total < qx > increases with invariant mass and 
bombarding energy in the same way as for hadron-hadron collisions EH H6 . 



13 



A similar behaviour of the average transverse momentum is found for central 
collisions of Pb + Pb at 10, 50 and 160 GeV/u in Fig. [IB], were we have plotted the 
total < q T > for both scenarios, i.e. without (dashed lines) and with (solid lines) 
'dropping' of meson masses. The channel decomposition is shown only for the 'free' 
mass scenario. Contrary to p + Be the contribution from 7i + tt~ annihilation dominates 
from 0.5 - 0.8 GeV at all bombarding energies, however, we do not find a pronounced 
difference in the total < qr > between both mass schemes. 



5.4 Rapidity distribution of lepton pairs 

The results of our calculations for the dilepton rapidity distribution in the nucleon- 
nucleon center-of-mass system are displayed in Figs. |17],|18| In analogy to the previous 
cases we study the systems p + Be at 10, 50, 450 GeV/u and Pb + Pb at 10, 50, 
160 GeV/u (b = 2 fm) integrating over the transverse momentum and the invariant 
mass range 0.4 < M < 0.7 GeV. As seen from Figs. |TT|,|r8| the shapes of the rapidity 



distributions are quite similar. An enhancement of about a factor 1.5-2 for Pb + Pb 
can be seen at midrapidity when employing in- medium meson masses (solid lines), 
while there is practically no difference at more forward or backward rapidities. Again 
this enhancement is most pronounced at 10 GeV/u for Pb + Pb and decreases with 
bombarding energy as expected from Fig. [I]. 



6 Direct photons 

Directly radiated thermal photons have been considered as an independent probe to 
study the hot and dense nuclear matter produced in ultrarelativistic nucleus-nucleus 



collisions |47], flq| . However, an experimental measurement of direct photons is a quite 
complicated task due to the background from hadronic decays. Only recently first 
upper limits for direct photon spectra have been reported by the WA80 collabora- 
tion for S + Au at 200 GeV/u. In the latter study photons from 7r° and r] Dalitz 
decays have been subtracted from the total photon signal; their spectra thus can be 
interpreted as upper bounds for the direct photon cross section. 

A first calculation of the direct photon radiation from a quark-gluon-plasma (QGP) 
was performed a couple of years ago in Ref. |50fl ; various hydrodynamical model 
calculations followed (cf. H9|, ^TJ and references therein), where the radiation from 
a QGP |52] has been compared to the radiation from a pure hadron gas scenario 



The comparison of the various models with the WA80 data P5fl, however, has 



demonstrated only the inapplicability of hadronic thermal models with high initial 
temperature. In this respect it should be quite useful to compare the WA80 upper 
limits with the results of a nonthermal model - such as the HSD transport approach 
- to find out possible conflicts with the hadronic scenario employed. 



6.1 Description of elementary channels 

In our analysis we take into account the following processes for photon production: 
ai — > 7T7, u — > 7T7, rj' — > fry or u^. The tc° and rj decays are already subtracted 
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experimentally and thus don't have to be taken into account in our calculations. 

The treatment of photon production in the HSD approach is quite similar to that 
for dileptons, however, using the branching ratios: 

~ 1.6 x 1(T 3 , = 0.085, = 0.3, = 0.03. (30) 

We discard baryon-baryon, meson-baryon and meson-meson bremsstrahlung in our 
present study since these channels were found to be of minor importance for dilepton 
production in case of S + Au at 200 GeV/u in Ref. p0[ . Furthermore, the soft- 



photon approximation employed in pDI is questionable at these energies according to 
the studies in Refs. @, . 

6.2 Comparison with experimental data for S + Au 

In Fig. [19| we show the result of our calculations for photon production in central 
S + Au collisions at 200 GeV/u in comparison with the experimental data |49|]. The 
computations were performed at b = 2 fm including the experimental rapidity cut 



2.1 < y < 2.9. As seen form Fig. [L^ the main contribution in our calculation comes 
from rf — > p/u 7 decays at low < 0.4 GeV and from uj — > Try for qx > 0.4; the solid 
line is the sum of all contributions which is still well below the upper limits of WA80. 

For the process a>i — > 777 we explore again both scenarios, i.e. without (dashed line) 
and with (dashed-dotted line) in-medium mass modification. In fact, the dropping of 
the a\ mass leads to a sizeable enhancement of a\ mesons in the reaction zone and 
thus to a significant enhancement of the photon spectra from the a\ as pointed out in 



Refs. |39], |51| . However, the relative contribution from the a± decay is still far below 
the 'background' from 7/ and uj decays. Thus, even in case of the 'dropping' meson 
masses we do not get in any conflict with the upper limits imposed by the WA80 data 



7 Summary 



In this work we have studied dilepton production in proton and heavy-ion induced 
reactions from 10 - 450 GeV or 10 - 160 GeV/u, respectively, on the basis of the 
covariant transport approach HSD [[| . We have incorporated the contributions of the 
Dalitz-decay of the rj,uj,T]',ai as well as 7r + 7r~ annihilation and the direct dilepton 
decay of the vector mesons p, uj, $ as well as K + K~ and 717? channels. It is found that 
for p + Be at 450 GeV the mesonic decays almost completely determine the dilepton 
yield, whereas in Pb + Au reactions the 7r + 7r~, K + K~ annihilation channels and 
np collisions contribute substantially. The experimental data taken by the CERES 



collaboration |TT], |43[ generally are underestimated by the calculations for invariant 
masses 0.35 GeV < M < 0.65 GeV when using 'free' form factors for the pion and 
p-meson in line with Refs. Jl9, 21], ^ . 

We have, furthermore, examined the effects of 'dropping' p, uj, $, K~ and a\ masses 
on the dilepton spectra, which generally leads to an improvement in the description of 
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the experimental data from the CERES [11], [|3| and HELIOS-3 collaborations (cf. also 
Refs. [[L9|, |2(| |29|). In order to allow for a clearer distinction between the different 
scanarios we have performed systematic studies for p + Be and central Pb + Pb 
collisions from 10 - 450 or 160 GeV/u, respectively, with respect to the transverse 
momentum distribution of the dilepton pair and their rapidity distribution. At all 
bombarding energies the in-medium effects are most pronounced at midrapidity, but 
the modifications in the transverse momentum spectra are only very moderate, such 
that the q± distributions do not qualify very much for quantifying the in-medium 
effects. The same holds for the average transverse momentum of the dilepton pair 
as a function of the invariant mass M. A more pronounced variation is obtained for 
the dilepton yield (integrated for M > 0.3 GeV) as a function of the charged particle 
multiplicity at midrapidity (cf. Fig. [7]), where the 'free' and 'dropping' meson mass 
scenarios differ significantly. This sensitivity should be explored in the next round of 
experiments at the SPS. 

Apart from the enhancement of the e + e~ yield for invariant masses 0.3 < M < 
0.7 GeV in case of 'dropping' meson masses, one should also see a decrease by factors 
4 - 6 in the latter scheme for 0.8 < M < 1 GeV in case of central Pb + Pb collisions 
(cf. Fig. [11]), provided that the experimental mass resolution is in the order of about 
10 MeV. The relative differences between the 'free' and 'dropping' mass schemes 
become larger when decreasing the bombarding energy from 160 GeV/u to 50 GeV/u 
and even to 10 GeV/u. According to Fig. |l] there might be a maximum sensitivity 
around 20 - 30 GeV/u, however, the relative change of the dilepton spectrum from 
10 - 50 GeV/u - apart from an overall increase by about a factor of 3 - is not very 
pronounced. 

We have, furthermore, explored if the hadronic reactions rates from the HSD 
transport approach might come in conflict with the upper limits for the direct photon 
spectrum in central collisions of S + Au from the WA80 collaboration P9| , because 
the photon multiplicity from a,\ decay might be dramatically enhanced in a phase with 
a partial restoration of chiral symmetry as suggested in Ref . [[51] . Our computations, 



however, including also the processes rf — > p/ui 7; 00 — > 7r7 do not show any indications 
for this, since the sum-spectra are still below the upper limits from WA80 by at least 
a factor of 3. In this respect we do not find any inconsistencies within the hadronic 
picture of the nucleus-nucleus collision. As pointed out in Ref. |56| , there might 
be problems due to the large energy densities reached in central Pb + Pb collisions 
(above 3 GeV/fm 3 ) at 160 GeV/u, where one generally believes that droplets of a 
quark-gluon-plasma (QGP) should be formed. 

We finally note, that apart from the multi-differential dilepton spectra - inves- 
tigated in this work - also dilepton angular anisotropies should provide additional 
information for disentangling the production channels experimentally. As analyzed in 
Ref. |33|] this additional observable, however, requires a large acceptance in the dilep- 
ton spectrometer, which is presently not fulfilled for the existing setups, but should 
be addressed in experiments with HADES. 
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Pb+Pb, central collision 
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Figure 1: The quantity F (Eq. (|25|)) for central Pb + Pb collisions as a function of 
the bombarding energy per nucleon for 4 different cuts in p m in- 
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Figure 2: The preliminary rapidity distribution of negative hadrons from NA49 |41 



(full squares) in comparison to the HSD results. The solid line corresponds to a 
calculation (at b = 2 fm) including the 'dropping' meson masses from Eqs. @, @, 
whereas the dotted line results from a calculation with bare meson masses. The open 
squares are obtained by reflecting the full squares at midrapidity. 
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Figure 3: Baryon density distribution (left column), momentum space (middle col- 
umn) and phase-space distribution (right column) for a 160 GeV/u Pb + Pb collision 
at b = fm for various times in fm/c. 
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Figure 4: The calculated dilepton spectra (full solid line) for p + Be at 450 GeV/u 
in comparison with the data from Ref. O]. The thin lines indicate the individual 
contributions from the different production channels including the CERES-acceptance 
and mass resolution; i.e. starting at low M: rj — > r ye + e~ (dashed line), uj — > 7r°e + e _ 
(dot-dashed line), rf — > •ye + e~ (long dashed line), a\ — > 7re + e _ (dot-dashed line); for 
M m 0.8 GeV: uj — > e + e _ (dot-dot-dashed line), p° — > e + e _ (dashed line), 7r + 7r~ — > 
e + e~ (dot-dashed line); for M ~ 1 GeV: $ — > e + e _ (dot-dashed line), 7rp — > e + e~ 
(dot-long dashed line), KK — > e + e~ (dashed line). 
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Figure 5: Dilepton invariant mass spectra for semicentral collisions of Pb + Au at 
160 GeV/u (full solid lines) in comparison to the preliminary data of the CERES 
collaboration ||43|| . The upper part shows the results of a calculation with bare meson 
masses whereas the lower part includes the 'dropping' meson masses (Eqs. (0),(0))- 
The assignment of the individual contributions is the same as in Fig. 4. 
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Figure 6: The charged particle multiplicity dn c h/dr] for Pb + Au at 160 GeV/u as a 
function of the impact parameter b for the free meson mass scheme (dashed line) and 
the dropping meson mass scenario (solid line). 
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Figure 7: The differential dilepton spectra for Pb + Au at 160 GeV/u integrated over 
the invariant mass region 0.3 < M < 1.0 GeV as a function of the charged particle 
multiplicity dn c h/drj without (open circles) and with (solid circles) the in-medium 
mass modification of the mesons. 
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Figure 8: The differential dilepton spectra for p + Be at 10, 50, and 450 GeV. The 
thick solid lines display the sum of all channels whereas the individual contributions 
are given in terms of the thinner lines. The assignment of the individual contributions 
is the same as in Fig. 4. The mass resolution employed is AM = 10 MeV. 



27 



E — r 



io° : 




— i — ■ — i — <— 

Pb+Pb, 10 GeV/u 

b=2 fm, free meson mass 




*-'-> I I i i fcfl — ^- . v 



0.6 



0.8 



1.0 



1.2 



E — r 



io° : 



£ io- 2 



3 10" r 



S^IO" 4 



10 



-5 ■ 



1 1 1 

Pb+Pb, 50 GeV/u 

b=2 fm, free meson mass 




co->7ie*e/'""*"~ 



I 

i 
i 



0.2 



0.4 



0.6 



0.8 



1.0 



1.2 



E — r 



io° : 



10 



1 ■ 



j „_„ + 



all 



1 1 1 1 — 

Pb+Pb, 160 GeV/u 

b=2 fm, free meson mass 



. 5 : \- 1 I >S 



0.2 



0.4 



0.6 



0.8 



1.0 



1.2 



M (GeV/c L ) 



Figure 9: The differential dilepton spectra for Pb + Pb at 10, 50, 160 GeV/u at 
b = 2 fm within the 'free' meson mass scenario. The assignment of the individual 
contributions is the same as in Fig. 4. The mass resolution employed is AM = 10 
MeV. 
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Figure 10: The differential dilepton spectra for Pb + Pb at 10, 50, 160 GeV/u at 
b = 2 fm within the in- medium meson mass scenario. The assignment of the individual 
contributions is the same as in Fig. 4. The mass resolution employed is AM = 10 
MeV. 
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Figure 11: The differential dilepton spectra for Pb + Pb at 10, 50, 160 GeV/u at 
b = 2 fm without (dashed lines) and with in-medium meson mass modification (solid 
lines) for a mass resolution AM = 10 MeV. 
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Figure 12: The channel decomposition for the average transverse momentum distri- 
bution in the invariant mass range 0.4 < M < 0.7 GeV for p + Be at 10, 50, and 
450 GeV; all contributions (solid line), uo — > 7r°e + e~ (dot-dot-dashed line), rj — >■ 7e + e~ 
(short dashed line), rj' — >■ 7e + e~ (long dashed line), p° — > e + e~ (dot-dashed line). 
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Figure 13: The transverse momentum distribution for the invariant mass range 0.4 < 
M < 0.7 GeV for Pb + Pb at 10, 50, and 160 GeV/u at b = 2 fm without (dashed 
lines) and with in- medium meson mass modification (solid lines). 
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Figure 14: The channel decomposition for the transverse momentum distribution for 
the invariant mass range 0.4 < M < 0.7 GeV for Pb + Pb at 160 GeV/u at b = 2 fm 
without (upper part) and with in-medium meson mass modification (lower part); all 
contributions (solid line), 7r + 7r~ — > e + e _ (upper dot-dashed line), uo — > 7r°e + e~ (short 
dot-dashed line), r\ — > 7e + e~ (dotted line), r/ — > 7e + e~ (dashed line), p° — > e + e _ 
(lower dot-dashed line). 
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Figure 15: The average transverse momentum < q T > (M) including the channel 
decomposition for p + Be at 10, 50 and 450 GeV. The assignment of the individual 
contributions is the same as in Fig. 4. 
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Figure 16: The average transverse momentum < > (M) including the channel 
decomposition for Pb + Pb at 10, 50 and 160 GeV/u at b = 2 fm without (dashed 
lines) and with in- medium meson mass modification (solid lines). The assignment of 
the individual contributions is the same as in Fig. 4. 
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Figure 17: The dilepton rapidity distribution for the invariant mass range 0.4 < M < 
0.7 GeV for p + Be at 10, 50, and 450 GeV. 



36 



10 1 , 



10° L 



10- 1 : 



10" 2 L 



10' 



— i 1 1 1 1 

Pb+Pb, 10 GeV/u 

b=2 fm, 0.4 < M < 0.7 GeV/c 2 , 




10 1 , 



10° L 



> 10 1 r 



10' 2 L 



10' 



12 
-| 1 1 1 1 1 — 

Pb+Pb, 50 GeV/u 

b=2fm, 0.4 < M < 0.7 GeV/c 2 



IP* 



tz 




free meson mass 

in-medium meson mass 

I i I i I i L 



io 1 F 







i 1 1 1 1 1 1 1 1 1 r 

Pb+Pb, 160 GeV/u 



b=2 fm, 0.4 < M < 0.7 GeV/c 2 




Figure 18: The dilepton rapidity distribution for the invariant mass range 0.4 < M < 
0.7 GeV for Pb + Pb at 10, 50 and 160 GeV/u at b = 2 fm without (dashed lines) 
and with in- medium meson mass modification (solid lines). 
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Figure 19: The differential photon multiplicity for central S + Au collisions at 
200 GeV/u in comparison to the upper limits from the WA80 collaboration [HUJ; 
all contributions (solid line). 
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